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We present a formalism to obtain equilibrium configurations of uni- 
formly rotating fluid in the second post-Newtonian approximation of gen- 
eral relativity. In our formalism, we need to solve 29 Poisson equations, 
but their source terms decrease rapidly enough at the external region of 
the matter(i.e., at worst 0(r -4 )). Hence these Poisson equations can be 
solved accurately as the boundary value problem using standard numer- 
ical methods. This formalism will be useful to obtain nonaxisymmetric 
uniformly rotating equilibrium configurations such as synchronized binary 
neutron stars just before merging and the Jacobi ellipsoid. 
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I. INTRODUCTION 



The last stage of coalescing binary neutron stars(BNS's) is one of the most promis- 
ing sources for kilometer size interferometric gravitational wave detectors, LIGO [l] and 
VIRGO 0. When the orbital separation of BNS's becomes ~ 700km as a result of the 
emission of gravitational waves, it is observed that the frequency of gravitational waves 
from them becomes ~ 10Hz. After then, the orbit of BNS's shrinks owing to the radiation 
reaction toward merging in a few minutes 0. In such a phase, BNS's are the strongly self- 
gravitating bound systems, and gravitational waves from them will have various general 
relativistic(GR) imformations. In particular, in the last few milliseconds before merging, 
BNS's are in a very strong GR gravitational field because the orbital separation is less 
than ten times of the Schwarzschild radius of the system. Thus, if we could detect the 
signal of gravitational waves radiated in the last few milliseconds, we would be able to 
observe directly the phenomena in the GR gravitational field. 

To interpret the implication of the signal of gravitational waves, we need to understand 
the theoretical mechanism of merging in detail. The little knowledge we have about the 
very last phase of BNS's is as follows: When the orbital separation of BNS's is < 10GM/ c 2 , 
where M is the total mass of BNS's, they move approximately in circular orbits because 
the timescale of the energy loss due to gravitational radiation tow is much longer than 
the orbital period P as 



where [l and d are the reduced mass and the separation of BNS's. Thus, BNS's adiabati- 
cally evolve radiating gravitational waves. However, when the orbital separation becomes 
6 — lOGM/c 2 , they cannot maintain the circular orbit because of instabilities due to the 
GR gravity @| or the tidal field ||. As a result of such instabilities, the circular orbit 
of BNS's changes into the plunging orbit to merge. This means that the nature of the 
signal of gravitational waves changes around the transition between the circular orbit 
and plunging one. Gravitational waves emitted at this transition region may bring us an 




(1.1) 
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important information about the structure of NS's because the location where the insta- 
bility occurs will depend on the equation of state(EOS) of NS sensitively |jU§. Thus, it 
is very important to investigate the location of the innermost stable circular orbit (ISCO) 
of BNS's. 

As mentioned above, the ISCO is determined not only by the GR effects, but also by 
the hydrodynamic one. We emphasize that the tidal effects depend strongly on the struc- 
ture of NS. Here, NS is a GR object because of its compactness, Gm/c 2 R ~ 0.2, where 
m and R are the mass and radius of NS. Thus, in order to know the location of the ISCO 
accurately, we need to solve the GR hydrodynamic equations in general. A strategy to 
search the ISCO in GR manner is as follows; since the timescale of the energy loss is much 
longer than the orbital period according to Eq. ( [Lip , we may suppose that the motion of 
BNS's is composed of the stationary part and the small radiation reaction part. From this 
physical point of view, we may consider that BNS's evolve quasi-stationally, and we can 
take the following procedure; first, neglecting the evolution due to gravitational radiation, 
equilibrium configurations are constructed, and then the radiation reaction is taken into 
account as a correction to the equilibrium configurations. The ISCO is determined from 
the point, where the dynamical instability for the equilibrium configurations occurs. It 
may be a grand challenge, however, to distinguish the stationary part from the nonsta- 
tionary one in general relativity. As Detweiler has pointed out |7]], a stationary solution 
of the Einstein equation with standing gravitational waves, which will be constructed by 
adding the incoming waves from infinity, may be a valuable approximation to physically 
realistic solutions. However, these solutions are not asymptotically flat because GWs 
contribute to the total energy of the system and the total energy of GWs inside a radius r 
grows linearly with r. The lack of asymptotic flatness forces us to consider only a bounded 
space and impose boundary conditions in the near zone. Careful consideration will be 
necessary to find out an appropriate boundary condition for describing the physically 
realistic system in the near zone. 

Recently, Wilson and his collaborators || proposed a simirelativistic approximation 
method in order to calculate the equilibrium configuration of BNS's just before merging. 
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In their method, they assume the line element as 

ds 2 = -(a 2 - PiP l )c 2 dt 2 + 2(3 i cdtdx i + ?p 4 dx 3 , (1.2) 

i.e., three metric 7^ is chosen as the conformal flat (i.e., 7^ = ip A 5ij), and solve only the 
constraint equations in the Einstein equation. In their approach, they claim that they 
ignore only the contribution of gravitational waves, but it is not correct at all; as shown in 
previous post-Newtonian(PN) analyses PHTTl, the tensor potential term exists in the three 



metric even if we ignore the radiation reaction of gravitational waves(i.e., V' _4 7ij 7^ <%). 
Since such a term appears from the second PN order in the PN approximation, the 
accuracy of their results is less than the 2PN order: In reality, from results by Cook et al. 
T2fl in which they obtain equilibrium configurations of the axisymmetric NS using both 



the Einstein equation and Wilson's method, we can see that some quantities obtained from 
Wilson's scheme, such as the lapse function, the three metric, the angular velocity, and so 
on, deviate from the exact solution by about 0((Gm/Rc 2 ) 2 ). This seems to indicate that 
their approach for the system of BNS's is valid only at the 1PN level from the PN point 
of view. Furthermore, the meaning of their approximation is obscure: It is not clear at all 
how to estimate errors due to such an approximation scheme and in which situation but 
the spherical symmetric system, the scheme based on the assumption of the conformal 
flatness is justified. 

In contrast with Wilson's method, the meaning of the PN approximation is fairly clear: 
In the PN approximation, the metric is formally expanded with respect to c" 1 assuming 
the slow motion and weak self-gravity of matter. If we will take into account the next PN 
order, the accuracy of approximate solutions will be improved. This means that we can 
estimate the order of magnitude of the error due to the ignorance of higher PN terms. 
Also, in the PN approximation, we can distinguish the radiation reaction terms, which 
begin at the 2.5PN order |I3| , from other terms in the metric. Thus, it is possible to 
construct the equilibrium configuration of BNS's without the radiation reaction terms in 
the 2PN approximation. 

We schematically describe two approaches in Tables 1(a) and 1(b). As mentioned 
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above, in close binary of NS's, it is important to take into account GR effects to orbital 
motion as well as to the internal structure of each NS. As for the orbital motion, there 
exist two parameters; one is the PN parameter v/c and the other is the mass ratio 77 of 
the reduced mass fi to the total mass M, and both parameters are less than unity. Thus, 
the physical quantities such as the orbital frequency are expanded with respect to them. 
In Table 1(a), we show schematically various levels of approximations in terms of v/c 
and rj. If all terms in a level are taken into account in the 2PN approximation, we mark 
P 2 N, while W means that all terms in the marked level are taken into account in Wilson's 
approach. From Table 1(a), we see that the 2PN approximation can include all corrections 
in 77 up to the 2PN order in contrast with Wilson's approach. On the other hand, Wilson's 
approach will hold completely in the test particle limit, i.e., at 0(r]°), whereas even in this 
limit the 2PN approximation is not valid at higher PN orders. As for the internal structure 
of each NS, there also exist two small parameters; one is the compactness Gm/c 2 R and 
the other is the deformation parameter from its spherical shape, such as an ellipticity e. 
In this case, the PN approximation becomes an expansion in terms of Gm/c 2 R. In Table 
1(b), we also show various levels of approximation in terms of these parameters. Although 
Wilson's approach is exact for spherical NS's, it is not valid in nonspherical cases even at 
the 2PN order. On the other hand, in the 2PN approximation, the spherical compact star 
cannot be obtained correctly in contrast with Wilson's approach. In this way, the 2PN 
approximation has a week point: Although it can take into account all effects up to the 
2PN order, it is inferior to Wilson's approach when we take a test-particle limit, 77 — > 0, 
or we describe an exactly spherical NS. However, as shown below, the error due to the 
ignorance of higher PN terms in those cases is not so large . 

To estimate the error due to the ignorance of the higher PN terms, let us compare the 
GR exact solutions with their PN approximations. First, we consider a small star of mass 
fx orbiting a Schwarzschild black hole of mass 3> /i. In this case, we may consider 
that the small star moves on the geodesic around the Schwarzschild black hole, and the 
orbital angular velocity becomes || 
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where r is the coordinate radius of the orbit in the harmonic coordinate. In the PN 
approximation, Eq. (|1.3|) becomes 



Comparing Eq.( |1.3| ) with Eq.( |1.4| ), it is found that the error size of the 2PN angular 
velocity is ~ 0.3% at f = 9GrribhC~ 2 , and ~ 1% at f = 6GmbhC~~ 2 . Thus, the 2PN 
approximation seems fairly good to describe the motion of relativistic binary stars just 
before coalescence. Next, we consider a spherical NS of a uniform density in order to 
investigate the applicability of the PN approximation for determination of the internal 
structure of NS's. In this model, the pressure, P, and the density, p =const., are related 



with each other 14 



P_ (1 - 2Gmr 2 /c 2 R 3 ) 1 / 2 - (1 - 2Gm/c 2 R) 1 ' 2 
~ 3(1 - 2Gm/c 2 Ry/ 2 - (1 - 2Gmr 2 /c 2 R 3 y/ 2 

1 Gm , r 2 s x G 2 m 2 , r 2 x Ghrf_ ,17 _ 19r| 3rj_ _ 8 
~2c 2 i?V R2j + C * R 2\ RJ>) + (£R3\g 8i? 2+ 8i? 4 8i? 6 ^ + ^ >' ^ ' 

where r s is the coordinate radius in the Schwarzschild coordinate and terms of order c~ 2 , 
c -4 and c~ 6 denote Newtonian, 1PN and 2PN terms respectively. In the second line in 
Eq. (|1.5|) , we expand the equation in power of Gm/c 2 R regarding it as a small quantity. 
In fig.l, we shows the error, 1 — P/P, in Newtonian, 1PN and 2PN function 
of r s for R = 5Gm/c 2 (solid lines) and 8Gm/c 2 (dotted lines), where P denotes the PN 
approximate pressure. It is found that the discrepancy in the Newtonian treatment is 
very large, while in the 2PN approximation the error is less than 10%. In this way, we 
can estimate rigidly the typical error size in the 2PN approximation. Furthermore, the 
accuracy is fairly good if the NS is not extremely compact; the 2PN approximation will 
be fairly accurate if the radius of NS is larger than ~10km. 

Thus, in the present paper, we develop a formalism to obtain equilibrium configura- 
tions of uniformly rotating fluid in the 2PN order as a first step. In section 2, we review 
the basic equations up to the 2PN order. In section 3, we rewrite the Poisson equation 
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for potential functions, which are described in section 2, into useful forms in which the 
source terms of the Poisson equations decrease rapidly enough(0(r~ 4 )). In section 4, 
we show a formulation to obtain numerically equilibrium solutions of uniformly rotating 
fluid in the 2PN approximation: Taking into account the formulation in the first PN 
approximation [J15|| , we further rewrite potentials defined in section 3 into a polynomial 
form in the angular velocity, Q. Then, we transform the integrated Euler equation into 
the polynomial form in Q 2 so that the convergence property in iteration procedures may 
be much improved. For the sake of analysis for numerical results, we describe the 2PN 
expression of the conserved quantities, such as the conserved mass, the ADM mass, the 
total energy and the total angular momentum in section 5. Section 6 is devoted to sum- 
mary. Throughout this paper, G and c denote the gravitational constant and the speed 
of light. Hereafter, we use units of G = 1. 



II. FORMULATION 

We write the line element in the following form; 

ds 2 = -(a 2 - pip^fdt 2 + 2(3 i cdtdx i + ip%jdx i dx j , (2.1) 

where we define det(7jj) = 1. To fix the gauge condition in the time coordinate, we use the 
maximal slice condition — 0, where K t l is the trace part of the extrinsic curvature, 
Kij. As the spatial gauge condition, we adopt the transverse gauge 7^ = in order 
to remove the gauge modes from 7^. In this case, up to the 2 PN approximation, each 
metric variable is expanded as |10 



1 U 1 
?2~ + ^ 

1 „ 1 /U 2 X I 



<A = 1 + ^ + -I(4)^ + 0(c- 6 ), (2.2) 



a = 1 - - 2 U + ? ( T + X) + - m a + 0(c- 7 ), (2.3) 

/3 J = ^(3)A + ^(5)A + 0(c- 7 ), (2.4) 

7« = S ij + ^h i] +0{c~ 5 ). (2.5) 

As for the energy-momentum tensor of the Einstein equation, we consider the perfect 
fluid as 



= (pc 2 + pe + P)u^u v + Pg^. (2.6) 

For simplicity, we assume that the matter obeys the polytropic equation of state(EOS); 

P = (r - l)pe = Kp r , (2.7) 

where T and K are the polytropic exponent and polytropic constant, respectively. Up to 
the 2PN order, the four velocity is expanded as fllBjliOl 



u° = 


1 + 


u = 






V 1 


It* = 






c . 




V* 


Ui = 






c 



1 /l 



1 /3 „.4 , 3 „.2rr , 1 rr2 



C 2v 2 ^ ~y 1 C 4V 8 ^ 1 2" ~ ' 2" ' -yj +°( c 

i + ^ 2 + ^) + Ul v " + ^ 2f/ + ^ 2 + - x )l + 



2 2 
c 4 V8 2 2 

5 2rr 1 

2 2 

where t>* = u l /u° and t> 2 = v l v % . Since we need w° up to 3PN order to obtain the 2PN 



+ v 1 ^ 4 + 7 -v 2 U + AU 2 -X + 4 (4 )^ + ( 3)/V') } + 0(c- 6 ), (2.8) 



equations of motion, we derive it here. Using Eq. (|2.8|) , we can calculate (aw ) 2 up to 3PN 
order as 



(cm ) 2 = 1 + tfj-^muj 



+ 2 (5) /V + (4(3)^^' + 4 ( 4)^ + yf/ 2 - 2X)v 2 + 8Uv 4 + v 6 } + 0(cT 7 ), (2 
where we use y % i = 5ij — c~^hij + 0(c~ 5 ). Thus, we obtain u° up to the 3PN order as 



= i + ^( V + u) + ±(§r,* + \v 2 u + V + (3) /V - x) 

+-{-( 6 )« + d(3)/3,-(3)/3,- + ZiyuV) + (5) ^y + 5 {3) (3 jV j U - 2UX 



c 2 ^2 / c 4 ^8 2 2 

1 

2 

^V-W^ ■ - W r ■ - 2 -J- ■ g - - lg 

Substituting PN expansions of metric and matter variables into the Einstein equation, 



+ (| (3 )/3y + 2 (4) V + 6f/ 2 - |x)t; 2 + yf/^ 4 + -^ 6 } + 0(cT 7 ). (2.10) 



and using the polytropic EOS, we find that the metric variables obey the following Poisson 
equations [|T0|1 ; 
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AC/ = -4vrp, 

AX = 4np(2v 2 + 2U + (3T - 2)e) , 
A (4 )V> = -27rp(V + e + -[/), 
A (3) f3i = lQnpv i - U h 



A (5)$ = 167Tp 



1 2 

(V + 2C/ + re) + {3) a - 4C/ !i ( (3) Aj + (3)/9i,j - *8mPk,k) 



_2 ( 4 ) ^ < + -(t7E/),i + (( 3 )A^i), 



A/ii,- 



-SijUAfiatU - SU^Uj + 5ijU t kU )k ) - 167r(/cwV - -<%/w 



2. * x _ , s 1 



-((3)A,j + (3)/3j,< - 3*<i(3)/3fc,fcJ - 2((X + 2 ( 4)^)^ - g-^-ApT + 2 (4) ^) 
A (6) a = 4tt P 2v 4 + 2v 2 (5C7 + re) + (3r - 2)eU + 4 (4) ^ + X + 4 (3) A? 
- - \uU tl Ui + Uj(2 {A) ^ - X) :l 

+ g ( 3 )^W (( 3 )^' + ( 3 )^ ~ 3^(3)^*.*) > 

where A is the flat Laplacian, and • denotes d/dt. 
Equations of motion for fluid are derived from 



v M rt = o. 



(2.11) 
(2.12) 
(2.13) 
(2.14) 

(2.15) 

(2.16) 



(2.17) 



(2.18) 



In this paper, we consider the uniformly rotating fluid around z-axis with the angular 
velocity Q, i.e., 



(2.19) 



where we choose fi J = (0, 0, Q) and is the completely anti-symmetric unit tensor. In 
this case, the following relations hold; 



(2.20) 



-dt dip 1 ^dt dip' ^dt dip' 

where Q, Qi and Qij are arbitrary scalars, vectors, and tensors, respectively. Then, 
Eq.( [2.18|) can be integrated as [17 

dP 



pc 2 + pe + P 



hlM U + C, 



(2.21) 



where C is a constant. For the polytropic EOS, Eq.( [2.2lD becomes 

YK 



In 



1 + 



c 2 (r - 1 



-p 



r-i 



or 



1 + 



TK 

c 2 (T-l) 



lnu° + C, 



u°exp(C). 



(2.22) 



(2.23) 



Using Eq. fl2.10|) , the 2PN approximation of Eq. (|2.22|) is written as 



H 2 H 



+ - 



"(6)° + o(3)A(3)A + 4 (3)/W^ 



?7 3 



+ (gjftuV + 2(4) V 



+ ^C/V + 2Uv 4 + \v & -UX- v 2 X + (5) /V + 



C, (2.24) 



where # = TKp^/iT - 1), v 2 = R 2 tt 2 and i? 2 = x 2 + y 2 . Note that Eq.(ggg) can be 



also obtained from the 2PN Euler equation like in the first PN case |18|, [15|. If we solve 
the coupled equations (2.11-17) and (|2.24|) , we can obtain equilibrium configurations of 
the non-axisymmetric uniformly rotating body. 



III. DERIVATION OF THE POISSON EQUATION OF COMPACT SOURCES 

FOR hij, ( 3) /3i and (5) /3j 

In section 2, we derive the Poisson equations for metric variables. However, the source 
terms in the Poisson equations for (3)/3j, (5) Pi, and fall off slowly as r — > 00 because these 
terms behave as 0(r~ 3 ) at r — > 00. These Poisson equations do not take convenient forms 
when we try to solve them as the boundary value problem in numerical calculation. Hence 
in the following, we rewrite them into other convenient forms in numerical calculation. 



As for h^, first of all, we split the equation into three parts as (10 
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A/ig = uiUtf - bijAU) - 3U ti Uj + SijUkUk = -4^$°, (3.1) 



A/iJf = -167r(pwV - \s ijP v 2 ), (3.2) 

2 

A^ G) = -((3) A j + (3)/3j,i - -5ij(3)$k,k) 

- 2((X + 2 (4) tP) tij - y. tJ A(X + 2 (4) ^)). (3.3) 
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The equation for h^p has a compact source, and also the source term of behaves as 



0(r" 



at r 



oo, so that Poisson equations for them are solved easily as the boundary 



value problem. On the other hand, the source term of h\ 3 - behaves as 0(r ) at r — > oo, 
so that it seems troublesome to solve the equation for it as the boundary value problem. 
In order to solve the equation for h\j as the boundary value problem, we had better 
rewrite the equation into useful forms. As shown in a previous paper 
integrated to give 

_ __ 

c 

1 



Eq.(|3J) is 



h 



2g- J(pv^x-y\d 3 y + 2—J(pv i y\x-y\d 3 y + 6,, 

I Plx - y ?d*y + J?— [ (pv 2 + W- P —) |x 
12 dx l dxJ J n y 1 y dx l dxi J v 2 / 1 



3 J x — v 



p|x - y\d y 

y\d 3 y 



(3.4) 



Using the relations 



p=-{po>\ i + 0(c-*) 



v l = 0, 



v l x l = 0, 



(3.5) 



Eq.( |3.4| ) is rewritten as 



+- 



■-x 

8 

1 



dx i J ^ x — y 



d 6 y + 



x-y 

d 



_d_ 

dx 



dxi 



(pv h y^jd 3 y 



t / (pv 2 + 3P - — ) a —d 6 y 4 

1 J v 2 / |x-y| y dxi 

2 f (pv 2 + 3P- pU/2 



3 3 



x- y 

2 ; |x-y| y 



y (pv 2 + 3P 



a- 



From Eq.( [3.6|) , it is found that h\f is written as 



h( (P = j( x %) p j + xj (s)Pi - QiP - Q% ') - Sij^mPk 



-x 



<9x 



■(^ {8) A-gg ) ) + 



+- 



(^ Q (')- Q m) + » (x'Q«-Qf>) -% Q m 



dxi 



2 
3< 



(3.6) 



(3.7) 
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where 



A( 3 )P = 


— 47T/w\ 




—Airpv^'x 3 , 


aq« = 


-4n(pv 2 + 3P 


aqP = 


-4vr(pt; 2 + 3P 



-p[/ x\ 



(3.8) 
(3.9) 
(3.10) 
(3.11) 



(G) 

Therefore, hh can be deduced from variables which satisfy the Poisson equations with 
compact sources. 

The source terms in the Poisson equations for ( 3 )/3, and (5)$ also fall off slowly. How- 
ever, if we rewrite them as fllTl 



)A = -4 (3 )Pi - -(a^-fc 

)A = -4( 5 )Pi - ^(2^(43-0 - T)i); 



where 



Aqi = —iTipx 1 , 
A (5) P = -4vrp 



1 2 

v\v 2 + 2U + Te) + (3) A + Uj^foj + (3) /?,-,; - ^y( 3 )/3*,k) 



(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 

then ( 3 )/3j and (5)$ can be obtained by solving the Poisson equations in which the fall-off 
of the source terms is fast enough, 0(r~ 5 ), for numerical calculation. Note that, using 
the relation ( 3 )P f = e^q^ and Eqs.( [2.20D , ( 3 )/3j and ( 5 )/3j may be written as 

( 3 )A = n{-4e izk q k + ]^(x l U^ - g^)} = fi (3 )A, (3.17) 

(5) A = fi{-4 (5) P + ^(2^(4)^ - (3-18) 

where 



A?]^ = -A-npiv 2 + e+ -U)x 



A (5) P = -4vrp 



"I * 2 * 

e izk x k (v 2 + 2U + Te) + (3) A + *7j((3)Aj + (3)A,i ~ o<%(3)A,fc) 



+ i(^),-^(( 3 )M fc ),, 



(3.19) 
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IV. DERIVATION OF BASIC EQUATIONS 



In this section, we derive the basic equation which has a suitable form to construct 
equilibrium configurations of uniformly rotating body in numerical calculation: Although 
equilibrium configurations can be formally obtained by solving Eq. fl2.24|) as well as metric 



potentials, U, X, ^ip, ^a, (3)A, (5)/?; and hij, they do not take convenient forms for 
numerical calculation. Thus, we here change Eq.( 2.24 ) into other forms appropriate to 
obtain numerically equilibrium configurations. 

In numerical calculation, the standard method to obtain equilibrium configurations is 
as follows iggiTJ; 



(1) We give a trial density configuration for p. 

(2) We solve the Poisson equations. 

(3) Using Eq.( [2.24| ), we give a new density configuration. 

These procedures are repeated until a sufficient convergence is achieved. Here, at (3), we 
need to specify unknown constants, Q and C. In standard numerical methods |19|j20|| , 
these are calculated during iteration fixing densities at two points; i.e., if we put p\ and 
p 2 at X\ and x 2 into Eq.( |2.24|) , they become two simultaneous equations for Q and C. 
Hence, we can calculate them. However, the procedure is not so simple in the PN case: 
VL is included in the source of the Poisson equations for the variables such as X, ^ip, 
( 6 )a, rji, (5)A, h\f\ QiJ\ and Q\ . Thus, if we use Eq. (|2.24j) as it is, equations for 



Vt and C become implicit equations for fl. As found in a previous paper |15| , in such 
a situation, the convergence to a solution is very slow. Therefore, we transform those 
equations into other forms in which the potentials as well as Eq.( |2.24j ) become explicit 
polynomial equations in Q. 

First of all, we define q 2 , q 2 i, <?4, q u , Qe and which satisfy 

Aq 2 = —AnpR 2 , (4.1) 

Aq 2i = -AnpR 2 x\ (4.2) 

Ag 4 = -4:7rpR 4 , (4.3) 

Aq u = -AnpU, (4.4) 
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Aq e = -Anpe, 
Aqij = —Ai\px % x 3 . 

Then, X, ^ijj, Q {1 \ Qf\ rji, ( 5) Pj, Q\J\ and h\p are written as 



(4.5) 
(4.6) 



v 


— ^q 2 \L zq u yoi ^)qe-, 


(A 7\ 


(4)V> 


= i(g 2 ^ 2 + g e + 7^q u ), 


(4.8) 




= q 2 n 2 + 3(r - l)q e - l -q u = q 2 tf + Q^, 


(4.9) 




= q 2i n 2 + Q$, 


(4.10) 


Vi 


= q 2i tt 2 + Voi, 


(4.11) 


(5) A 


= Cizk&kQ 2 + (5)-Poi) 


(4.12) 






(4.13) 




= 4TL 2 {e iz kej Z iqki — -<%g 2 ), 


(4.14) 



where Q$ , Voi and < b )P Qi satisfy 



AQ 



(J) 

Oi 



-4tt(3P - -pU)x l = -47rp(3(r - l)e - ^[/ 



Ar] 0i = -A-K P (e + hj)x\ 
A (5 )Poi = -4?rp e i2fc a; fc (2f/ + re) + (3) ^ + C/j( (3 )A,j + (3)&\i - -<%3)4,fc) 



2 
-( 

3 



(4.15) 
(4.16) 



(4.17) 



(GO 

Note that (5)$ and /i^- are the cubic and quadratic equations in il, respectively, as 

1 



(5) A = ^ 



-4 (5) Poi + ^{^(?e + j^u)^ - 

(5) A (A) fi + (5) A (B) ^ 3 , 



(4.18) 



(G) 



9 



_d_ 

dx % 



If g 

2la^ 



T £jzkQk,ip X J * £izkQk,ip £izkQkj,tp £jzkQki,ip^J &ijX ^kzlQl 



-A— 



-x 



8 19a;' 



9 



^ £kziQi,tp — £kziqij,tp 



) + dxi (' 



^ ekziqi^ — tkziQi 



%,ip 



(4.19) 
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Finally, we write as 

(6) a = (6)ao + (6)tt2^ 2 - 2g 4 f2 4 , 
where ( 6 )« and ( 6 )« 2 satisfy 
A (6) a = Anp 



(4.20) 



(3T-2)eU- (3T-A)q e + 3q u 

3 
2 



■(/4f + - -t^tr, + u^fa + (3r + i)g e ) 



-47r5 (ao) , 



(4.21) 



A (6) a 2 = %-npR 2 (hU + Te + 2( 3 )^) - (Ae izk e jz iq M - -<%g 2 + h[f^U# + Sq 2 ,iUi 



1 - / . . 2 



= -4vr5 (a2) . 



Using the above quantities, Eq. (|2.24j ) is rewritten as 

tt2 tt3 p6 

# - ^ + — A = a + Bn 2 + Dn A + ^-n 6 + c, 

2c 2 3c 4 6c 4 ' 

where 



(4.22) 



(4.23) 



A—U + ^{2q u + (3r - 2)g e ) + ^{- (6 )«o - ^ + C/(2g u + (3r - 2)g e ) }, 
i? 2 1 11 

5 = y + ^( 2R2U + 2 <?2 + (3)4) + ^{-(6)«2 + 2(3)A(3)A + 4(3)4^ 
+ (3r - l)q e R 2 + \q u R 2 + ^U 2 R 2 + 2q 2 U + (5) /^ + \{h^ + 

D = ^ + ^{ 2 <? 4 + (3)4^ + ^2^ 2 + 2f/^ 4 + (bj/JW + I(/ig> + AR 2 q RR )}. (4.24) 

Note that in the above, we use the following relations which hold for arbitrary vector Qi 
and symmetric tensor Qij, 

Qip yQx ~t~ ^cQyi 

Q<p<p 11 Q XX 2<XyQ ' X y ~\- X Qyyi 

R 2 Qrr = x 2 Qxx + ZxyQxy + y 2 Q yy - (4.25) 

We also note that source terms of Poisson equations for variables which appear in A, B 
and D do not depend on Q explicitly. Thus, Eq. Q4.23Q takes the desired form for numerical 
calculation. 
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In this formalism, we need to solve 29 Poisson equations for U, q x , q y , q z , (5)Pq x , (5)Po y , 



q yz , (6)Q;o, (6)Q!2 and g 4 . In Table 2, we show the list of the Poisson equations to be solved. 
In Table 3, we also summarize what variables are needed to calculate the metric variables 
U, X, ( 4 )^, (e) a, (3) A, (5)A, f47\ ^ and ^f } - Note th at we do not need ( 5 )Po z , ?7o 2 , 
and g 22 because they do not appear in any equation. Also, we do not have to solve the 
Poisson equations for hJ£p and q yy because they can be calculated from h z u ^ = —h^—h^p 
and q yy = q 2 - q xx - 

In order to derive U, qi, q 2 , ?2i, <?4, 9e and g^-, we do not need any other potential 
because only matter variables appear in the source terms of their Poisson equations. On 
the other hand, for q u , Qcii\ Voi and h\f\ we need the Newtonian potential U, and for 
(5)-fbi) (6) a o and (6)«2, we need the Newtonian as well as PN potentials. Thus, U, qi, q 2 , 
Q2i, q4, qe and qij must be solved first, and then q u , QoPj Voi, hf^\ (5)Poi and (6) a 2 should 
be solved. ( 6 )tt must be solved after we obtain q u because its Poisson equation involves 
q u in the source term. In Table 2, we also list potentials which are included in the source 
terms of the Poisson equations for other potentials. 

The configuration which we are most interested in and would like to obtain is the 
equilibrium state for BNS's of equal mass. Hence, we show the boundary condition at 
r — > oo for this problem. When we consider equilibrium configurations for BNS's where 
the center of mass for each NS is on the rr-axis, boundary conditions for potentials at 
r — > oo become 



Vox, Vo y , Qox , Qo y , Qoz\ 92, Q2x, Q2 y , fe, 9«, 9e, hf x \ hf y \ h x u z \ hf y \ hf z \ q xx , q xy , q. 



XZ1 




(4.26) 




(5)-M)y 



^Jsp X dV + ^fspydV + 0(r-% 
^JSy P) -dV + ^jS y ^ydV + 0(r-% 



(4.27) 



16 



n y 

^y = ^ 



n -\ px 2 {e + \u)dV + 0^)., 
fpy 2 (e + ^u)dV + 0(r- 4 ), 



(4.28) 



o {I) 

Vox 
Qoy 
Qoz 



% J px 2 {3(T - l)e - -U)dV + 0{r 
^ J py*(3(r-l)e-±u)dV + 0(r- 
^ | pz 2 (3(T-l)e-^u)dV + 0(r- 



g 2 , = ^ JpR 2 x 2 dV + 0(r-*), 
n y r 

q 2y = -j pRYdV + 0(r^), 

q 2z = ^JpR 2 z 2 dV + 0(r- 4 ), (4.29) 



IL yy 



l -js^dV + 0{r-% 

J SfJyzdV + 0(? 



3n y n 



3n x n y 



3n x n z 



S^xydV + 0(r 
SifxzdV + 0{ 



-5^ 



(4.30) 



</,, - - / px 2 dV + 0(r ••). 
3n x n 2 



(6)«o = J /^dK + 0(r- 3 ), 



where dV = d 3 a;, and 



^3/ = „a px y dV + <J{r 



3n y n z 



py 2 z 2 dV + 0(r" 5 ), (4.31) 



(6) a 2 = i y S^dV + 0(r" 3 ), (4.32) 



We note that at r -»• oo, -> 0(r- 5 ), -> 0(r- 6 ), -> 0(r~ 4 ) and 

0(r -4 ), so that all the above integrals are well defined. 



(4.33) 



(a 2 ) 



V. CONSERVED QUANTITIES 

In this section, we show the conserved quantities in the 2PN approximation because 
they will be useful to investigate the stability property of equilibrium solutions obtained 
in numerical calculations. 
(1) Conserved mass 
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J p*d 3 x, 



where 



(5.1) 



0„/,6 



= pau ip 
= P 



1 + -^ 2 + 3U) + + ^ 2 C7 + ±V + 6 (4) V + (3)/W) + 0(c" 6 ) 



c 2 ^2 / c 4 ^8 2 

Equation ( |5.2| ) may be written as 



P* = P 



1 /l 



1/3, 13 



45, 



1 + + 3f/ ) + ?(r 4 + y v2u + ~4 u2 + 3Ue + ^ iVi ) + 0(c-6) 



(2)ADM mass |2T|,[T0 



Madm = -g- y A^rf 3 x = y p ADM d 3 x, 



• (5-2) 



(5.3) 



(5.4) 



where 



Padm — P 



1 2 1 

+ 2(3)/^ + (3) ((3) A j + (3)Pj,i ~ 3^(3)^*,*)} + °( c " 6 ) 



, (5.5) 



or 



Padm — P 



1 + i(v a + e + -£/) + 1(V + 9t; 2 t/ + iW + 5Efe + ^[/ 2 + || (3) /V 

+0(0-' 



(5.6) 



(3)Total energy, which is calculated from Madm ~ in the third PN order [10 

E = J p E d 3 x, 



(5.7) 



where 
Pe = P 



\v 2 + e- \u) + i (V + \v*U + Tvh + 2Ue - \u* + \ {3) frV 
1 rll . .4/^ , 47 rA , , fiTwr , 41 r2 5 



I 16 + V re + + ' 4(4) ^ + 6rw + ¥ + 2 (3)/w " x 

- V + 2T m f3 i v i e + 5e ( 4)^ + 5*7 (3) /V - y^(4)^ + ^ 

+-hijV l V J + -( 3 )/3i(3)A 



+ 



167rp 



(zhijU^j + (3)A',i((3)A,i + (3)^?,i - g^j(3)A,fc)^| + 0(c 



(5.8) 
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(4) Total linear and angular momenta: In the case — 0, these are calculated from [21 

1 



Pi 



lim f Kijirv'dS 

iff r— >oo 



— lim <p ib^Kjjni dS 



1 



57T 



v 107T 7 



(5.9) 



where Jj = (pc 2 + pe + P)au°Ui. Up to the 2PN order, the second term in the last line of 
Eq. (pT9|) becomes 

1 



167T 
1 

167T 
1 

IQtt 



hjk,i(3)Pj,kd x, 



d 3 x, 



Jim 4 h jk ,i(3)Pjn dS = 0, 



(5.10) 



where we use — > 0(r x ) and (3)/3j — > 0{r 2 ) at r — > oo, and the gauge condition 



hjk k — 0. Thus, in the 2PN approximation, Pj becomes 



where 
Pi = P 



P% = J Pid 3 x, 



+ ^{ vi ( y2 + Te + 6U ) + (3)^} + ^{ V + (5) A + (3)A(v 2 + 6U + Te 



(5.11) 



+v i (2 (3) /3 i v i + 10 ( 4) ^ + 6IW + — U 2 + Yev 2 + IWv 2 + v A -X)\ + 0(c 



The total angular momentum J becomes 



J = J Pipd^x, 



(5.13) 



(5.12) 



where p v = -yp x + xp y . 



VI. SUMMARY 



It is generally expected that there exists no Killing vector in the spacetime of coalescing 
BNS's because such a spacetime is filled with gravitational radiation which propagates to 
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null infinity. However, we may consider coalescing BNS's as the almost stationary object 
from physical point of view as described in section 1. Motivated by this idea, in this 
paper, we have developed a formalism to obtain equilibrium configurations of uniformly 
rotating fluid up to the 2PN order using the PN approximation. The concept of being 
"almost" stationary becomes clear in the framework of the PN approximation and, in 
particular, the stationary rotating objects can exist exactly at the 2PN order, since the 
energy loss due to the gravitational radiation does occur from the 2.5PN order. There 
appear, at the 2PN order, tensor potentials which were completely ignored in Wilson's 
approach ||. It should be noted that these tensor potentials play an important role at 
the 2PN order: This is because they appear in the equations to determine equilibrium 
configurations as shown in previous sections and they also contribute to the total energy 
and angular momentum of systems. This means that if we performed the stability analysis 
ignoring the tensor potentials, we might reach an incorrect conclusion. 

In our formalism, we extract terms depending on the angular velocity Q from the 
integrated Euler equation and Poisson equations for potentials, and rewrite the integrated 
Euler equation as an explicit equation in Q. This reduction will improve the convergence 
in numerical iteration procedure. As a result, the number of Poisson equations we need to 
solve in each step of iteration reaches 29. However, source terms of the Poisson equations 
decrease rapidly enough, at worst 0(r -4 ), in the region far from the source, so that we 
can solve accurately these equations as the boundary value problem like in the case of the 
first PN calculations ||15|| . Thus, the present formalism will be useful to obtain equilibrium 



configurations for synchronized BNS's or the Jacobi ellipsoid. These configurations will 
be obtained in future work. 
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Figure Captions 

Fig. 1 Error of the pressure in the post-Newtonian approximation for the GR compact 
star of uniform density as a function of the normalized areal radius(r/i?). Solid 
and dotted lines show the case R = bGm/c 2 and 8Gm/c 2 , where R and m are the 
circumference raduis and the mass of star, respectively. 
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Table 1 (a) 

Various levels of approximation in terms of PN expansions(t> 2 /c 2 ) and mass ratio(r? = 
/jl/M; /j, = reduced mass, M = total mass). We mark P 2 N if all terms in that level are 
taken into account in the 2PN approximation, while W is marked if Wison's approach 
takes into account all terms in that level. The mark — means that the relevant term does 
not exist and the levels taken into account by neither approaches are blank. We neglect 
secular effects due to gravitational radiation reaction in Tables 1(a) and (b). It should 
be noted that, at 0{jf), Wilson's approach produces exact GR solutions, but it is not 
justified at the 2PN order even at 0{rj r ). 



PN \ Tj 


V° 


v 1 


v 2 


o(v 3 ) 


N 


P 2 N, W 








1PN 


P 2 N, W 


P 2 N, W 






2PN 


P 2 N, W 


P 2 N 


P 2 N 




>3PN 


W 









Table 1 (b) 

Various levels of approximation in terms of PN expansions(Gm/c 2 i?) and ellipticity of 
a NS(e). The meanings of P 2 N and W are the same as those in Table 1(a). Wilson's 
approach produces exact GR solutions in the case of the completely spherical star. 



PN \ e 


e = 


e ^ 


N 


P 2 N, W 


P 2 N, W 


1PN 


P 2 N, W 


P 2 N, W 


2PN 


P 2 N, W 


P 2 N 


>3PN 


W 





Table 2 



List of potentials to be solved(column 1), Poisson equations for them(column 2), and 
other potential variables which appear in the source term of the Poisson equation(column 
3). Note that i and j run x,y,z. Also, note that we do not have to solve r] 0z , (s)Po z , q yy , 
q zz and hf z \ 



Pot. 


Eq. 


Needed pots. 




Pot. 


Eq. 


Needed pots. 


U 


(2.11) 


None 




Qij 


(4.6) 


None 


Q% 


(3.14) 


None 




{I) 


(4.15) 


U 




(4.1) 


None 




Voi 


(4.16) 


u 


q2i 


(4.2) 


None 






(4.17) 


U, qi 


Qa 


(4.3) 


None 




(6)«0 


(4.21) 


Una h {U) {I) 


Qu 


(4.4) 


U 




(6)«2 


(4.22) 


U, q 2 , Qi, Q2i, Qij 




(4.5) 


None 






(3.1) 


U 



Table 3 

Variables to be solved in order to obtain the original metric variables. 



Metric 


Variables to be solved 


see Eq. 


U 


U 


(2.11) 


(3) Pi 


Qi-, U 


(3.17) 


X 


12, Qu, Qe 


(4.7) 




<?2, Qu, Qe 


(4.8) 




(5)-f()i) VOi, Qu, Qe 


(4.18) 


(5)A-' B) 


Q2i, Q2 


(4.18) 


(6)« 


(6)«0, (6)«2, QA 


(4.20) 


h iU) 


h iU) 


(3.1) 


h (S) 


Qij,Q2 


(4.14) 


h (A) 


Qoi, Qu, Qe 


(4.19) 


h {B) 


Qij, Q2, Q2i, Qi 


(4.19) 



